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ABSTRACT 

9C J1503+4528 is a very young CSS radio galaxy, with an age of order 10 4 years. This 
source is therefore an ideal laboratory for the study of the intrinsic host galaxy/IGM proper- 
ties, interactions between the radio source and surrounding ISM, links between star formation 
and AGN activity and the radio source triggering mechanism. Here we present the results 
of a spectroscopic analysis of this source, considering each of these aspects of radio source 
physics. 

We find that shock ionization by the young radio source is important in the central re- 
gions of the galaxy on scales similar to that of the radio source itself, whilst evidence for an 
AGN ionization cone is observed at greater distances. Line and continuum features require the 
presence of a young stellar population, the best-fit model for which implies an age of 5 x 10 6 
years, significantly older than the radio source. 

Most interestingly, the relative sizes of radio source and extended emission line region 
suggest that both AGN and radio source are triggered at approximately the same time. 

If both the triggering of the radio source activity and the formation of the young stellar 
population had the same underlying cause, this source provides a sequence for the events 
surrounding the triggering process. We propose that the AGN activity in 9C J1503+4528 was 
caused by a relatively minor interaction, and that a super-massive black hole powering the 
radio jets must have been in place before the AGN was triggered. 

Key words: galaxies: active - galaxies: evolution - galaxies: ISM - galaxies: stellar content 
- radio continuum: galaxies - galaxies: individual: 9C J1503+4528 



1 INTRODUCTION 

The triggering of powerful radio sources is still poorly understood. 
It is generally thought that some mechanism must channel signifi- 
cant quantities of gas into the central regions of the host to fuel the 
AGN, wherein a rapidly spinning black hole may power the radio 
jets (e.g Blandford and Znajek, 1977). 

Two scenarios have been proposed. The first is the infall of 
gas which has cooled in the halo of a galaxy, or a cooling flow in 
a cluster (Bremer, Fabian & Crawford 1997); this proposition has 
received much interest recently because it may have a thermostat 
feedback effect to prevent runaway cooling in clusters (Briiggen 
and Kaiser 2002, Fabian et al. 2003). However it is most likely that 
this mechanism in fact causes a renewal of activity in an existing 
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radio source (Nipoti and Binney 2005), rather than being related to 
the original triggering. 

The second is associated with the long-standing idea that a 
merger between two galaxies, at least one of which is gas rich, is a 
precursor to AGN activity (e.g. Heckman et al 1986). Wilson and 
Colbert (1995) suggested that the coalescence of black holes fol- 
lowing a major merger creates the necessary rapidly- spinning black 
hole, which suggests links between radio sources and other merger 
systems such as ULIRGs, and implies that evidence for young stel- 
lar populations and disturbed morphologies should be present at 
some level in observations of young radio sources. 

The merger scenario is supported by an increasing amount of 
anecdotal evidence. For instance, there are several cases known 
where the host galaxy of a powerful radio source appears to be 
undergoing a major merger (e.g., Johnston et al. 2005). Similarly, 
there is increasing evidence that powerful high-redshift radio galax- 
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ies lie in forming clusters (see e.g. Venemans et al. 2005); such en- 
vironments would be rich in major mergers. This evidence in turn 
takes us back to the first scenario above; massive cluster galax- 
ies may host a spinning black hole as a relic of their early history, 
which is generally quiescent but which becomes active when fu- 
elled by cooling gas. It is clearly evident that radio sources play a 
significant part in hierarchical galaxy and structure formation (e.g. 
Croton et al 2006), so the precise ingredients for radio source for- 
mation are a cause of renewed interest. 

In any case, our understanding of this early phase of radio 
galaxy evolution is hampered by the difficulty of compiling uni- 
form samples of radio sources in the first few millennia of their 
growth, and indeed the impossibility of pinpointing galaxies im- 
mediately prior to this stage. However, it is is now widely accepted 
that Compact Steep Spectrum (CSS) and some Gigahertz Peaked 
Spectrum (GPS) radio sources represent the earliest stages of radio 
source evolution, being the progenitors of the larger FRII source 
population (e.g. Fanti 1995, Readhead 1996, Snellen et al 2000, 
Murgia 2003). The study of such sources in the stages immediately 
after triggering may be a valuable way of elucidating the events that 
have occurred beforehand. 

In an effort to secure a uniform sample of young sources, we 
have turned to the 9C survey (Waldram et al 2003). This was car- 
ried out at 15 GHz with the Ryle Telescope to identify foreground 
sources for the Very Small Array (VSA) microwave background 
experiment. Over the past several years, extensive followup of the 
radio properties of this population has been carried out (Bolton et 
al 2003, 2004, 2006; Bolton et al submitted), both to understand 
the effects on microwave background experiments better and to in- 
vestigate the sources themselves; the high survey frequency has 
proved a good way of selecting significant numbers of CSS and 
GPS sources, as well as sources peaking in the ~ 10 GHz range. 
We are now undertaking optical followup of the CSS and GPS 
sources, and in this paper we present a detailed spectroscopic study 
of one of these compact sources, 9C J1503+4528, a classical dou- 
ble FRII radio galaxy at a redshift of z ~ 0.5 with an angular size 
of 0.5". 

Details of the observations are presented in section 2, analysis 
of the emission line strengths and ionization state in section 3, gas 
kinematics in section 4, emission morphology in section 5 and a 
discussion of the implications of our results in section 6. We give a 
summary of our conclusions in section 7. 

Values for the cosmological parameters of Qo = 0.27, Qa = 
0.73 and H = 65km s^Mpc" 1 are assumed throughout. 



case of 9C J1503+4528 implies that the source has an age of or- 
der ~ 10 4 years. This places the source very close to its time of 
triggering, since the typical lifetime of an FRII source is ~ 10 7 to 
10 8 years (e.g. Kaiser, Dennett-Thorpe and Alexander 1997). Fi- 
nally, the source redshift (z ~ 0.5) makes many useful diagnostic 
lines and the 4000-A break accessible to moderate-resolution spec- 
troscopy in a wavelength range region relatively clear of night sky 
emission. With these facts in mind, the young CSS radio source 
9C J1503+4528 is therefore a particularly appropriate object for 
the study of the intrinsic host galaxy/IGM properties, interactions 
between the radio source and surrounding ISM, links between star 
formation and AGN activity and the underlying radio source trig- 
gering mechanism. 



2.1 VLA imaging 

The 22-GHz data obtained by Bolton at el. (2004) were re-analysed 
to obtain the map shown in Fig. Q Both the central component 
and the SE compact component are unresolved. The most obvi- 
ous interpretation of the radio morphology is that the bright cen- 
tral component is the core, with the SE compact component be- 
ing a hotspot, and the SE and NW extensions being lobes. Given 
the overall steep spectrum, it is implausible that the source is be- 
ing viewed at an angle close to the line of sight; in such a case, a 
strongly boosted flat-spectrum jet component would be expected 
to dominate in this high-frequency map. We thus conclude that 
the source has a double-lobed structure with an overall extent of 
~ 3kpc. 



2.2 Optical and infrared imaging 

An image using the Isaac Newton telescope Wide Field Camera r- 
band filter was taken on the night of 2005 April 11. This filter has 
the advantage of not including the strong [Oil] and [Om] emission 
lines at the source redshift. Five 600-s images were taken, offset at 
the centre and corners of a 20" pattern, and stacked. The resulting 
image is shown in Fig.Q 

A Jf-band image was also obtained, on the night of 2006 
March 5 using the UFTI instrument at UKIRT. Nine 60-s jittered 
observations were taken and reduced using the standard ORAC-DR 
pipeline; the resulting image is also displayed in Fig.Q The see- 
ing, measured from an observation of UKIRT standard star P177-D 
taken immediately after the target observation, was 0.6 arcseconds. 



2 THE OBSERVATIONS 

The observed source, 9C 1503+4528 (see Fig.[TJ, was selected as 
our first target for detailed study for several reasons. 

First, the source has a rest-frame 1.4 GHz luminosity of 1.33 
xlO 25 W Hz sr -1 , well above the FRI/FRII break. Second, the 
source is resolved, and clearly shows a classical double structure in 
the radio map presented in Bolton et al. (2004). The angular size 
of the source is 0.5", corresponding to a projected linear size of 
~ 3 kpc. It is currently accepted that GPS and CSS sources are at 
a young stage in the evolution of the radio source (see e.g. Snellen 
et al 2000). The picture in which sources are young is strongly sup- 
ported by observations of the hotspot advance speed in many of 
the most compact sources. These observations find hotspot advance 
speeds in the range 0. 1-0. 4c (e.g. Polatidis and Conway 2003; Tay- 
lor et al 2000; Owsianik, Conway & Polatidis 1998), which in the 



2.3 Optical Spectroscopy 

The spectroscopic observations were carried out on 2004 July 18, 
using the Keck LRIS spectrograph (Oke et al 1995). Seeing was 
~ 0.8", and a 1" slit was used. The 600/7500 grating was used, 
providing a wavelength range of ~ 5200 - 7800A, a spatial scale 
of 0.215 arcsec/pixel, a spectral scale of 1.26A/pixel, and a spec- 
tral resolution of ~ 4.4A. For an object lying at a redshift of 
z ~ 0.5, the corresponding rest-frame spectral range is approx- 
imately 3500 — 5200A, thereby including emission lines from 
[Oll]3727A up to [Om]5007A. Slit PAs of +45° and -45° were 
chosen, aligning the slit either parallel or perpendicular to the radio 
source axis. Long-slit spectra were obtained with total integration 
times of 1800s at each slit PA. The integrations were split into 900-s 
exposures to aid the removal of cosmic rays. 

Standard packages within the NOAO IRAF reduction software 
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Figure 1. The compact radio source 9C J1503+4528. Left: UFTI K— band (top) and WFC r— band (bottom) images of the host galaxy, an R ss 19 elliptical. 
The two spectroscopic slit PAs, parallel and perpendicular to the radio source axis are indicated on each diagram by the solid and dashed lines respectively. 
Centre: VLA 22GHz radio contours. Right: VLA and OVRO radio spectrum for this source, from Bolton et al (2004). 



were used to reduce the raw data. Corrections were made for over- 
scan bias subtraction, and the data were then flat-fielded using in- 
ternal lamps. The flat fields were taken immediately after each ob- 
servation of the target, to minimize the effects of fringe drift in the 
CCD. The two-dimensional frames were wavelength calibrated and 
corrected for distortions in the spatial direction using wavelength- 
calibration arc exposures. Taking care not to subtract any of the 
extended line emission, the sky background was removed. Obser- 
vations of the smooth-spectrum standard star BD +28 4211 were 
used to provide accurate flux calibration of the final spectra and to 
remove the telluric A- and B-band absorption features. The central 
regions of the parallel and perpendicular spectra effectively probe 
the same region of the host galaxy. These two spectra were com- 
bined in order to produce a higher signal-to-noise spectrum, for 
use in analysing the emission from the central regions of the radio 
galaxy. 

The resulting parallel and perpendicular 2-d spectra are dis- 
played in Fig. [2] In the parallel spectrum, the presence of narrow 
emission, extended over a distance of ~ 5" (approximately 10 
times larger than the radio source itself) can clearly be seen, in addi- 
tion to a broader component lying close to the continuum centroid. 
The [Oil] emission is particularly strong compared to 03, and the 
continuum emission short of the 4000A break is also substantially 
more luminous than would be expected for an old stellar popula- 
tion alone. The perpendicular spectrum displays similar properties, 
with the exception that extended narrow emission is not observed. 
A 1-d spectrum extracted from the central 2" of the combined per- 
pendicular and parallel spectra is displayed in Fig. [5] 

Determining the precise systemic redshift of 9C J 1503+4528 
was complicated by the inherent structure in the emission lines: the 
redshift implied by the peak of the line emission on the contin- 
uum centroid differs substantially from that of the narrow extended 
component. However, the various absorption features in the spectra 



(Ca H & K, high-order Balmer lines) all suggest a redshift identical 
to that determined from the extended narrow component, giving a 
value of z = 0.521. 



3 EMISSION LINE FLUXES AND THE IONIZATION 
STATE OF THE GAS 

One-dimensional spectra were extracted from the parallel, perpen- 
dicular and combined 2-d frames. The extractions were carried out 
over 2" (~ 12kpc) for the parallel, perpendicular and combined 
frames, and additionally over 5" (~ 31kpc) for the parallel and 
perpendicular frames. The lines were then fitted by gaussians; the 
resulting fluxes and equivalent widths are presented in Appendix 
A. 

Spectra were also extracted in 0.5" steps across the slit, in or- 
der to investigate the changing ionization state of the emission line 
gas with position via ionization state diagnostic diagrams (Bald- 
win, Phillips & Terlevich 1981). For an effective analysis, model 
predictions should be well separated and results should be consis- 
tent over several different diagrams/line combinations. Line pairs 
should ideally be close in wavelength and of the same species, to 
limit the effects of reddening and variations in gas composition. 
The most appropriate lines available from our optical spectra are 
[On]3727A, 034363A, H/3 and 035007 A. 

Figure [4] displays the diagnostic diagram for the line pairs 
[Oil] 3727 A/H/3 and 035OO7A/H/3, extracted from the parallel 
spectra in 0.5" steps, together with various theoretical model 
tracks. These include: (i) the shock models of Dopita & Suther- 
land (1996), which also incorporate the effects of a precursor ion- 
ization field, (ii) the mixed medium (matter bounded vs. ionization 
bounded) photoionization tracks of Binette et al (1996), plus (iii) 
simple AGN photoionization tracks (Groves et al 2004a, 2004b), 
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Figure 2. Sections of particular interest (rest frame 3730-4070A and 5190- 
5490A) from the 2-d optical spectra of 9C1503+4528, with the slit aligned 
parallel (left) and perpendicular (right) to the radio source axis. Wavelength 
increases from the bottom to the top of the diagram, as indicated by the 
direction of the arrow, and the position of the 4000A break is denoted by 
the dashed line. The extracted regions are approximately 10" wide. The five 
prominent emission lines are (bottom to top): [OlI]3727A, [NeIII]3869A, 
H/3, 034959A and 035OO7A. 



also incorporating the effects of dust. The photoionization models 
displayed are for solar metallicity and a density of n = 1000cm -3 . 
For reference, the photoionization tracks display little variation for 
lower values of n, while the n = 10000cm -3 track is shifted to 
smaller values of log([OlI]/H/3) by ~ 0.2dex. The 2Z Q metallicity 
tracks are generally similar to those of solar metallicity, whereas 
lower Z tracks are shifted to smaller values of both log([Oll]/H/3) 
and log(03/H/3). 
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Figure 3. Spectrum of 9C J1503+4528 obtained from the central 2" of the 
combined parallel and perpendicular data. This spectrum provides the high- 
est signal-to-noise ratio for our data, and provides a good illustration of the 
relative emission line strengths and the fairly blue colour of the continuum. 



This diagram clearly shows the changing ionization state rel- 
ative to distance from the nucleus. The spectra extracted from the 
central regions of the host galaxy lie close to the predictions of the 
shock ionization models (and also the far end of the mixed medium 
ionization track). At intermediate distances (~ 1", i.e. sampling 
emission from just outside the radio source), the shock plus precur- 
sor photoionization models are most suited. Beyond 1-1.5" (i.e. the 
extended narrow emission region), the data are well described by 
simple AGN photoionization. High values of the ionization param- 
eter, U, are the most plausible, and a dusty medium is favoured. For 
the photoionized gas, the ionization state at larger radii suggests ei- 
ther a lower spectral index (a ~ —2.0 rather than a ~ —1.2), 
higher gas density or lower metallicity. The mixed medium mod- 
els of Binette et al only provide a good explanation for a few of 
the data points, and it seems unlikely that this system can be well 
described using that particular model. The observed line ratios are 
not likely to be affected by internal reddening, as the Balmer line 
ratios in our spectra are in good agreement with the predictions 
of case B recombination. Balmer line absorption could potentially 
have some influence on the exact positioning of the data points in 
this figure; while we believe we have correctly accounted for the 
errors on our measurement of H-/3, additional shifts of up to 0.2dex 
towards smaller line ratios cannot be completely ruled out. How- 
ever, this would not have any effect on our overall interpretation of 
the data. 

The picture provided by these results is as follows. The UV 
continuum emission from the AGN ionizes the gas at all locations 
within the ionization cone, decreasing in intensity at larger dis- 
tances. The innermost regions of the galaxy, of order the same size 
as the fledgling radio source, have an additional ionization con- 
tribution from shocks. The UV radiation field produced by these 
shocks adds another ionization component to the emission from 
gas at intermediate distances. 

The relative ionization states of the parallel and perpendicular 
spectra can be compared in Fig. [5] This figure displays the same 
tracks as in Fig.|4] for spectra extracted from various regions of our 
2-d spectra. These include the central 2" regions of the combined, 
parallel and perpendicular spectra. The results for the central 1-2" 
should not be greatly different, as both spectra sample essentially 
the same region. However, the impact of AGN and shock ioniza- 
tion will be greatly lessened at larger distances, for material which 
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Figure 4. Ionization mechanism diagnostic plot for 9C1503+4528, using 
the line ratios [OlI]3727/H,3 and [OlII]5007/H/3. Data points represent 
spectra extracted in 0.5" steps, and are labelled in terms of their distance 
(in arcseconds) from the continuum centroid. The model tracks include: (1) 
shock ionization (Dopita & Sutherland 1996), with and without a precursor 
photoionization region (light & dark shading respectively), (2) the matter- 
bounded photoionization track (dashed dark blue track, Binette, Wilson & 
Storchi-Bergmann 1996), and (3) simple AGN photoionization tracks (light 
blue tracks), including the effects of dust (solid black tracks; Groves, Do- 
pita & Sutherland 2004a,b). The photoionization models are for n=1000 
cm -3 , solar metallicity and various values of the ionization index a (from 
1.2 < a < 2.0). A brief description of the effects of changes in metallicity 
and/or gas density is given in the main text; we refer the reader to Groves et 
al (2004a,b) for full details. 

lies off the radio source axis and well outside the ionization cone 
of the AGN. Also displayed are the observed emission line ratios 
for spectra extracted at distances from the continuum centroid of 
+ 1.5" to +2.5" and -1.5" to -2.5", for both the parallel and 
perpendicular spectra. The large errors on the data points for the 
extended regions of the perpendicular spectra reflect the very weak 
line emission present at these distances. In the case of the —1.5" 
to —2.5" data point, shocks are not a plausible option, and the data 
should lie outside the ionization cone of the AGN. Instead, the po- 
sition of this data point suggests that ionization by young, hot stars 
may provide a better explanation (see e.g. Kewley, Geller & Jansen 
(2004) and Rola, Terlevich & Terlevich (1997) for appropriate lo- 
cations on the [Oll]/H/3 vs. 03/H/3 diagnostic plot for Hll regions; 
such emission line ratios would generally be expected to lie towards 
the lower left quadrant of Fig. [5}- A combination of ionization by 
young stars and the AGN could also be relevant for the —1.75" and 
+2.25" data points on Fig.|4] 

Figure |6| shows the diagnostic diagram for the 
035OO7A/034363A and [Oll]3727 A/035007 A line ratios, to- 
gether with the observed data from our parallel spectra in spatial 
steps of 0.5". Note that the mixed medium model track, which 
bent towards the data points at low values of Ajy// 1 in Fig. [4] 
tends towards higher values of 035007 A/034363A on Fig. [6] as 
opposed to the lower values of this line ratio observed in the data. 
Although the error bars are large for the 035007 A/034363A line 
ratio, we nevertheless have a certain amount of consistency with 
the previous diagnostic diagrams: once again the central regions 
of the galaxy favour the predictions of shocks, while more distant 

1 The solid-angle ratio of matter-bounded to ionization-bounded material 
in the mixed-medium model of Binette et al 1996. 
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Figure 5. Ionization mechanism diagnostic plots for 9C1503+4528, us- 
ing the line ratios [OlI]3727/H/3 and [OlII]5007/H/3. Data points represent 
spectra extracted from various regions of the spectra taken perpendicular 
and parallel to the radio axis, and include: (i) the central 2" of the paral- 
lel, perpendicular & combined (best signal-to-noise) spectra, (ii) 1" wide 
sections of the parallel and perpendicular spectra lying ±2" from the con- 
tinuum centroid. The model tracks are the same as those in Fig.RTi. 
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Figure 6. Ionization mechanism diagnostic plots for 9C1503+4528, using 
the line ratios [OlI]3727/[OlII]5007 and [OlII]5007/[OlII]4363. Data points 
represent spectra extracted in 0.5" steps, and are labelled in terms of their 
distance from the continuum centroid. The model tracks are the same as 
those in Fig.|4]i, and include: (1) shock ionization (Dopita & Sutherland 
1996), with and without a precursor photoionization region (light & dark 
shading respectively), (2) the matter-bounded photoionization track (dashed 
dark blue track, Binette, Wilson & Storchi-Bergmann 1996), and (3) Sim- 
ple AGN photoionization tracks (light blue, various styles), including the 
effects of dust (solid black tracks, Groves, Dopita & Sutherland 2004a,b). 

locations lie closer to the dusty AGN photoionization tracks. 
The relatively low 035007 A/034363A ratio is indicative of high 
temperatures and/or densities in the emitting gas, of the order of 
several 10 4 K in the low density limit (e.g. Osterbrock 1989), fully 
consistent with the presence of shocks in these regions. 



4 GAS KINEMATICS 

We have also investigated the changing gas kinematics, via an anal- 
ysis of the profiles of the [Oll]3727A and 03 (4959+5007 A) lines. 
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Figure 7. Example line profile fits for the 03(4959+5007)1 emission lines. 
The top figure is for the spectrum extracted at a distance of 2.25" to the SE 
of the continuum centroid, and the bottom figure for a region lying at a dis- 
tance of 0.75" in the same direction. The data were modelled with several 
components (dotted lines): a continuum level, and 1-3 separate Gaussian 
components for each line. For 03, the relative strengths and positions of the 
4959A and 5007 A lines were fixed. The combination of gaussians with the 
lowest reduced \ 2 determined the maximum number of velocity compo- 
nents which could be fitted to the data. 
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Figure 8. Gas kinematics for 9C J1503+4528. The changing line fluxes 
(top), widths (centre) and offsets (bottom) with position are displayed 
for the [Oil] (left) and 03 (right) emission lines, extracted and fitted in 
0.5"steps from the parallel spectrum. Up to three emission line compo- 
nents are modelled at each position. Negative offsets are towards the NW, 
and positive offsets towards the SE. 



Two-dimensional regions around each of these emission lines were 
extracted, and from these a sequence of one-dimensional spectra 
were extracted along the slit direction, stepped every 0.5". 

After allowing for continuum subtraction, the data were fitted 
by up to three successive gaussian components, which were only 
accepted if their FWHM was larger than the instrumental resolu- 
tion and the signal-to-noise ratio was greater than approximately 
five (excepting the outermost spectra, for which a fainter, single 
gaussian component could be detected at a lower signal-to-noise.) 
The combination of gaussians with the lowest reduced \ 2 were 
then taken to be the maximum number of velocity components 
which could be fitted to the data (we also fixed the relative posi- 
tions and strengths of the 4959A and 5007A lines). Example fits 
are displayed in Fig. [7] 

The integrated line flux, the velocity relative to that at the cen- 
tre of the galaxy, and the emission line FWHM were determined 
for each gaussian fit. We calculated the FWHM by subtracting in 
quadrature the instrumental FWHM, as determined from unblended 
sky lines. We then calculated errors for these three parameters, al- 
lowing for the fact that a range of possible fits are equally accept- 
able at low S/N. This approach allows us to search for high velocity 
components in the emission line gas, or other structures incompat- 
ible with a fit to a single velocity component. Given the observed 
emission line profile of this source, such a step was a necessity, and 
broader emission line components are clearly present in the central 
parts of the extended emission line region. 



The results of our analysis of the kinematics of the parallel 
spectrum are presented in Fig. [8] 

An unresolved emission component is present in the outer re- 
gions of the extended emission line region (hereafter EELR), in 
both [Oil] and 03, at wavelengths consistent with the emitting ma- 
terial lying at rest relative to the host galaxy. This narrow com- 
ponent can be traced across the galaxy in 03. The [Oil] emission 
line is of substantially higher luminosity, particularly in the central 
regions of the galaxy. Because of this, the narrow component of 
[Oil] can only be observed in the outer regions, and is lost in the 
centre in favour of one or more broader components. The relative 
strengths of the narrow components of [Oil] and 03 are consistent 
with simple AGN photoionization, as demonstrated by our diag- 
nostic diagrams in the previous section. 

The broader emission line components extend to approxi- 
mately 1-1.5" from the continuum centroid. Note that the [Oil] 
emission is substantially stronger than the 03 emission; such a low 
ionization state in the emitting material is suggestive of shocks as- 
sociated with the radio source, in which case the emitting material 
would display larger line widths than more settled material. Usually 
a single broad component is sufficient, except in the central regions 
where two broad components are preferred. The broad components 
on the SE side of the galaxy (see the 2-d spectrum presented in 
Fig. |2j usually lie at shorter wavelengths, with a velocity shift of 
a few hundred km s~\ whilst the emission from the opposite side 
of the galaxy is shifted to longer wavelengths by a similar amount. 
Of these two distinct components, the blue-shifted emission line 
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Figure 9. Gas kinematics for 9C J1503+4528. The changing line fluxes 
(top), widths (centre) and offsets (bottom) with position are displayed for 
the [Oil] (left) and 03 (right) emission lines, extracted and fitted in 0.5" 
steps from the perpendicular spectrum. Up to three emission line compo- 
nents are modelled at each position. Negative offsets are towards the SW, 
and positive offsets towards the NE. 



component is generally more luminous, and has a larger line width 
(400-600 km s~\ cf 200-400km s" 1 for the red-shifted emission). 

We also present the results of a similar analysis of the spec- 
trum aligned perpendicular to the radio source axis (Fig. |9}. As 
the line emission in the perpendicular spectrum is much less exten- 
sive, larger extracted spectra of 1" in width have been modelled to 
boost the signal-to-noise. The first main feature of note is the lack 
of an extended narrow component. This is strong evidence that the 
source of the ionizing photons for the extended narrow component 
is the ionization cone of the AGN, which does not intersect with 
the perpendicular spectrum except in its centre. 

In the central regions of the perpendicular spectrum, broad 
emission components are observed. These are very similar to those 
in the parallel spectrum, with the exception that there is no clear 
tendency for the blue- or red-shifted emission to preferentially lie 
on one side of the galaxy or the other. This suggests that outflows 
along the radio axis are a plausible option for explaining the gas 
kinematics, as an alternative to simple rotation. The greater lumi- 
nosity and extent of the material to the SE side of the galaxy, to- 
gether with the blueshifted broader component, suggests that the 
SE radio lobe would be orientated towards the observer; given the 
observed radio asymmetries (section 2.1) and quiescent extended 
emission line component at larger radii the outflow scenario is in- 
deed highly plausible. 



5 CONTINUUM EMISSION 

The continuum emission of 9C 1503+4528 is remarkably blue. To- 
gether with the clear presence of B aimer absorption lines and the 
young age of the radio source this suggests a significant young stel- 
lar population. If the triggering of radio source activity is connected 
with the star formation, the properties of any young stellar popu- 
lation formed as part of that process can shed light on the exact 
mechanisms and timescales involved. We have therefore carried out 
continuum modelling of our spectra, with the aim of placing limits 
on both the mass and age of any young stellar population present in 
this galaxy. 

For the continuum modelling, we combined both the parallel 
and perpendicular 2-d frames and extracted a 1-d spectrum from 
the central 2" (Fig.|3j. 

The first stage was to remove the contribution of nebular con- 
tinuum emission (Dickson et al 1995). Although subject to stel- 
lar absorption features, the Balmer line strengths are in general 
consistent with the predictions of case B recombination, and in- 
ternal reddening is unlikely to be a major issue. The H/3 flux can be 
well modelled, and was used to determine the contribution from the 
various nebular continuum emission processes (free-free emission, 
bound-free emission, two-photon emission and the Balmer forest 
pseudo-continuum), assuming an electron temperature of 15000K. 
(Checks with alternative temperatures, allowing for the error on our 
evaluation of the H/3 flux, do not produce significantly different re- 
sults.) The corrected spectrum was binned in blocks of ~ 30A, 
excluding any regions where emission lines were present. 

To fit the remaining stellar continuum we used a combina- 
tion of young and old stellar population SEDs from the models of 
Bruzual & Chariot (2003). A Chabrier initial mass function and 
solar metallicity were selected, and the stellar populations were as- 
sumed to form in an instantaneous burst. We consider stellar pop- 
ulation ages ranging from < 10 5 years up to 9 x 10 9 years (the 
expected age of the old stellar population in a z ~ 0.5 galaxy, as- 
suming a formation redshift of z ~ 10 — 20). Reduced \ 2 fitting 
was used to determine the relative proportions of old/young stars 
for all combinations of SED models. 

Figure [Tol displays the results of our best-fit model. This con- 
sists of an old (6 Gyr) stellar population, plus a 5 x 10 6 year young 
stellar population accounting for 0.065% of the total stellar mass, 
and roughly 70% of the total continuum flux at a wavelength of 
~ 3600A. Our modelling implies that the total stellar mass in the 
central region of the galaxy probed by our spectra is approximately 
10 OMq, whereas the overall stellar mass of the galaxy (based on 
our if— band modelling in section 6) is closer to IO^Mq. The 
largest residuals lie near the Ca absorption feature at ~ 3970A, 
which is subject to infilling by line emission at 3967A due to He + 
[Nem]. 

A plot of best-fit reduced \ 2 as a function of young stellar 
population (hereafter YSP) age is displayed in Fig. ll II A contribu- 
tion from a young stellar population is very much favoured over an 
old stellar population alone. Fits obtained using different ages for 
the majority old stellar population (ranging from 5 to 9 Gyr in age) 
do not lead to any significant variation in the resulting best-fit YSP. 

We must also consider the potential presence of a power-law 
ionizing continuum associated with the AGN. If the presence of 
a power-law component is simply added to our modelling, allow- 
ing its flux and shape to vary freely, the resulting best-fit model 
generally prefers an unrealistically strong power-law flux, and it is 
therefore necessary to place constraints on the power law contri- 
bution in order to make a meaningful assessment of its effects. To 
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Figure 10. Continuum fits to the central 2" of the combined spectrum (Fig.|3j. The top frame (a) displays the observed spectrum (dotted line), the best-fit 
model (dark solid) and the two component (old and young) SEDs which make up the best-fit model. Also displayed in frame (a) are the model residuals for 
each of the bins used in the fitting procedure. For clarity, we also display the observed data (frame (b), centre left), the best-fit model (frame (c), lower left), 
the old stellar population component (a 6 Gyr old stellar population, frame (d), centre right) and the young stellar population component (a 5 million year old 
stellar population, frame (e), lower right) as separate spectra. 



determine the balance between YSP age/mass and power-law con- 
tribution, we repeated our fitting procedure using fixed power-law 
contributions ranging from 0% to 50% of the total flux at a wave- 
length of A ~ 3600A. The spectral index of the power law (defined 
as Fx <x X a ) is constrained to be one of three values, a — — 1, 
—2 or —3. The fitting procedure was then repeated for all combi- 
nations of stellar populations and power-law parameters. The re- 
sulting best-fit YSP parameters (age, mass, reduced x 2 ) for each of 



the different power law contributions considered are displayed in 
Appendix B. 

With the exception of the fits with the very highest power-law 
contribution (> 25%), the best-fit YSP mass and age remain rela- 
tively unchanged. The similar strengths of the narrow line emission 
components either side of the galaxy centroid in the parallel spec- 
tra (suggesting that this source is observed with the ionization cone 
aligned in a direction nearly perpendicular to the line of sight) also 
imply that a two-component stellar population provides a more than 
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Figure 11. Reduced x 2 of best-fit model vs. young stellar population com- 
ponent age for the continuum fitting of 9C J1503+4528. The best-fit model 
includes a YSP component with an age of 5 X 10 6 years, accounting for 
0.065% of the total stellar mass of the host galaxy, and roughly 70% of 
the observed continuum flux at a wavelength of ~ 3600A. The age of the 
majority old stellar population in the best-fit model is approximately 6Gyr. 

adequate explanation of the observed continuum emission, and that 
no power-law is required. 



6 IMPLICATIONS FOR SOURCE TRIGGERING 

6.1 Morphological properties of the host galaxy 

Given the spectroscopic evidence for a young stellar population, 
we examined the galaxy morphology for evidence for a merger or 
other interaction which may have caused the starburst. Our calcu- 
lated 0.06% 5 Myr YSP would be expected to contribute only 0.02 
magnitudes to the A-band flux, and so the A"-band image will be 
dominated by the old stellar populations of the host galaxy and any 
merging companion or tidal tails. 

We fitted the galaxy morphology using the methods detailed 
in Inskip et al. (2005), which can be briefly summarised as follows. 
The point spread function for the UFTI field is determined from un- 
saturated stars within the mosaic field; we extract their 2-d profile, 
normalise to unit flux and take the average profile. Sersic profile 
(Sersic 1968) galaxy models were then convolved with the point 
spread function and fitted to the surface profile of the galaxy, using 
available least squares minimisation IDL routines 2 . Free parame- 
ters are the galaxy flux, centroid, effective radius, Sersic index and 
fractional nuclear point source contribution. Galaxy ellipticity was 
also allowed to vary, but in the case of this object a circular model 
with e = was preferred. 

From the A-band data, our resulting best-fit model, residu- 
als and 1-d radial profile are displayed in Fig. 1121 We find that the 
host galaxy appears to have a simple de Vaucouleurs elliptical mor- 
phology, with a preferred Sersic index of 4.0 ± 0.2, and a best-fit 
(X 2 = 1-03) effective radius of r cff = 1.70" ± 0.15" (equivalent 
to 11.5 ± l.Okpc at z = 0.521). We find no evidence for any nu- 
clear point source contribution at A" band (preferred point source 
percentage < 1%). 

2 MP2DFUNFIT. PRO, part of Craig Markwardt's MPFIT non-linear least 
squares curve fitting package available via http://astrog.physics.wisc.edu/ 
~craigm/idl/fitting.html. 




Figure 12. Host galaxy morphology fits for 9C J1503+4528. The K-baad 
galaxy image is presented in frame (a), with the model galaxy given in 
frame (b) and the residuals after subtracting the best fitting model from the 
data displayed in frame (c). Frame (d) illustrates a 1-d cut through the data 
(dotted line and binned points) with the best-fit model (solid line) overlaid. 

The r-band data are too noisy to make an accurate fit to the 
data; the effective radius is constrained only to lie between 1.3" ± 
0.5" and 1.6" ± 0.5". This is, however, consistent with the A- 
band fit. Most importantly, we find no significant residual emission 
which would indicate the presence of a close companion galaxy or 
tidal tails. 

Any galaxy interaction/merger activity must therefore either 
have been very minor in nature, or have occurred sufficiently long 
ago that the old stellar population of the host galaxy has returned to 
an undisturbed morphology. In the latter case, the interaction must 
have taken place several dynamical timescales before the present 
observations. We estimate the dynamical timescale Td yn as r e ff/er, 
where a is the velocity dispersion of the galaxy, which in turn we 
estimate by taking the FWHM of the narrow component of the 
emission lines (~ 200 km s _1 ). This gives Td yn ~ 60 Myr, so 
any previous interaction must have occurred at a time very much 
earlier that the formation of the YSP. We thus conclude that the 
YSP and present AGN activity were caused by, at most, a minor 
interaction. 

6.2 Evidence from the extended emission-line region 

As with many other radio galaxies the properties of the extended 
emission line region of 9C J 1503+4528 vary substantially with spa- 
tial position. The changes in ionization state confirm our expecta- 
tions that shocks are important in the vicinity of the young radio 
source, and that AGN photoionization dominates elsewhere. The 
kinematic properties of the emission line gas confirm this picture, 
with the broader emission lines highlighting the disturbed kinemat- 
ics of the central regions, in contrast to the kinematically quiescent 
material at larger radii. 

The projected physical size of the EELR is ~ 34 kpc. Emis- 
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sion line regions of this scale are by no means unusual, and are 
often substantially larger in the case of FRII radio sources (for 
example, ~50-100 kpc for extended 3C and 6C sources, Best, 
Rottgering & Longair 2000; Inskip et al 2002). However, there 
are significant contrasts between this EELR and those of more ex- 
tended radio galaxies. 

First, the EELR of 9C J1503+4528 has a very sharp cutoff, 
rather than the more gradual decline in emission line luminosity 
observed in the cases of other sources (e.g. Best, Rottgering & Lon- 
gair 2000; Inskip et al 2002). Second, and of great significance, is 
the fact that the EELR is a factor of ten larger than the size of the 
radio source. Unlike the case of the 3C and 6C sources, where the 
EELR is significantly affected by the passage of the radio source, 
here we have a case where the EELR could very well correspond 
to the limiting distance reached by the ionization front of an AGN, 
illuminating the region ahead of the radio source before it leaves 
the centre of the host galaxy. 

This factor of ~ 10 difference between radio source size and 
EELR size has a significant implication. In combination with the 
typical 0.1 c source expansion speed, we find it highly plausible 
that the radio jets and AGN were triggered simultaneously, some 
~ 10 4 years prior to our observations. 

6.3 Spatial distribution of star formation 

The material lying at larger angles to the radio source axis (i.e. 
the regions probed by the perpendicular spectrum) is also of in- 
terest. The line emission along this axis is far less extensive (as ex- 
pected for the regions lying on the edge of, or outside, the ionization 
cone). Its ionization state suggests a lower ionization parameter U 
(i.e. fewer and/or less energetic ionizing photons), and the weak 
line emission in these regions could be better explained as being 
due to ionization by young OB stars (HII regions). Our continuum 
modelling in the preceding section finds clear evidence for a young 
stellar population with an age of roughly 5 million years, but does 
not allow us to pinpoint its location. However, the varying strength 
of the B aimer absorption features in the parallel and perpendicular 
spectra CFig. 1131 can give us some idea of how widely distributed 
this young stellar population might be. 

This plot displays the H/3 line profile for the parallel and per- 
pendicular spectra, using increasing extraction widths of l"to 5", 
beyond which point increased aperture widths only lead to a deteri- 
oration of the signal-to-noise. In the perpendicular extracted spec- 
tra, Balmer absorption can be clearly observed, and becomes in- 
creasingly apparent as emission from regions at a greater distance 
from the AGN is included (i.e. for increased aperture widths). The 
line emission in the parallel spectra is roughly 50-100% more lumi- 
nous (see Appendix A), and the evidence for the same trend in the 
strength of the absorption component is less clear-cut. Whilst the 
line emission is strongest in the central regions of the host galaxy, 
the distribution of the line absorption must therefore be less well 
peaked, otherwise we would not observe a strengthening of the 
absorption component in wider apertures. This suggests that any 
young stellar population responsible for the Balmer absorption fea- 
tures is unlikely to lie predominantly within the central few pc. The 
pronounced appearance of the Balmer absorption in the perpendic- 
ular spectra suggests that the underlying young stellar population 
is not linked to the location of the growing radio source but rather, 
more widely distributed throughout the galaxy. Indeed, given the 
relative ages of YSP and radio source (5 x 10 6 years and ~ 10 4 
years respectively), a radio-source induced star formation scenario 
would seem unlikely. 
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Figure 13. H/3 line profiles for the parallel (left) and perpendicular (right) 
spectra, extracted using increasing aperture widths of (from top): 1", 2", 
3", 4" and 5". Balmer absorption is most clear in the perpendicular spectra, 
and generally becomes more pronounced with the inclusion of emission at 
greater distances from the AGN (i.e. the larger extraction apertures). 



It is interesting to contrast the derived properties of the YSP 
detected in the spectra of 9C J1503+4528 with those found for other 
sources. Recent work on the GPS radio source PKS1345+12 (Ro- 
driguez et al, in prep) finds evidence for a very similar young stellar 
population to that of 9C J1503+4528: an age of roughly 5 million 
years (c.f. ~ 10 3 years for the radio source itself) and a mass of 
~ 10 6 Mq. Imaging observations locate the YSPs in several super 
star clusters lying in the outer regions of the host galaxy, rather than 
in a nuclear starburst ring such as that observed for Cygnus A (e.g. 
Jackson, Tadhunter & Sparks 1998). A further point of interest is 
the finding that the hosts of high luminosity AGN observed as part 
of the Sloan Digital Sky Survey (Kauffmann et al 2003) typically 
have younger mean stellar ages than normal elliptical galaxies ex- 
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isting at the same epoch. Additionally, the young stars are not pref- 
erentially located near the galaxy nuclei, but rather are spread out 
over a much wider scale. These properties are strikingly similar to 
those observed in 9C J1503+4528. 



6.4 Gas kinematics 

The observed changes in gas kinematics and dominant ionization 
mechanism in different regions of the EELR are not unexpected, 
given the expected location of radio source shocks, but do provide 
two key results: the clumpy ISM/IGM responsible for the EELRs 
observed around powerful radio sources does appear to be in place 
prior to the expansion of the radio source, and the kinematic prop- 
erties of the extended regions are relatively undisturbed. Any dis- 
turbance of the ISM/IGM due to an earlier interaction/merger has 
since settled down. 



6.5 Spectroscopic evidence of a prior merger 

If the triggering of radio source activity in this galaxy was due to a 
merger or interaction, the properties of any young stellar population 
formed as part of that process can shed light on the exact mecha- 
nisms and time scales involved. Our continuum modelling of these 
spectra confirms the presence of a young stellar population, whose 
formation occurred several million years in advance of the onset 
of AGN activity. The mass of this YSP component is a low pro- 
portion of the total stellar mass of the host galaxy. Together with 
the the quiescent kinematics in the outer regions of the EELR, and 
the morphological evidence described above, this suggests that if a 
merger/interaction was responsible for both the eventual triggering 
of the radio source activity and the formation of the young stellar 
population, it is most likely to have been fairly minor in nature. 



7 CONCLUSIONS 

9C J1503+4528 provides a detailed insight into many of the events 
associated with the triggering of a powerful radio source. Our ob- 
servations have shown that: 

• The relative sizes of EELR and radio source (the former being 
« 10 times larger), together with typical compact radio source ex- 
pansion velocities of 0. 1 c, suggest that both the luminous quasar 
nucleus responsible for photoionizing the outer regions of the 
EELR, and the radio source itself, were triggered simultaneously. 

• The EELR ionization state is consistent with shocks being im- 
portant on scales comparable to the radio source, with AGN pho- 
toionization within the ionization cones being the dominant mech- 
anism further from the AGN, beyond the extent of the radio lobes. 

• The gas kinematics confirm that weak shocks could be im- 
portant close to the radio source, and that the outer regions of the 
EELR are relatively undisturbed. The kinematics suggest that mo- 
tion along the radio axis may be taking place, which may indicate 
the presence of outflows. 

• The blue continuum emission and Balmer absorption features 
suggest that recent star formation has occurred, but their spatial 
distribution suggests that it is not caused by the passage of the radio 
source. 

• Modelling of the observed continuum emission gives a best-fit 
age for the young stellar population of 5 x 10 6 years, and a mass 
of roughly 0.06% of the total stellar mass of the host galaxy. 



• Imaging observations of the host galaxy reveal a simple ellipti- 
cal morphology; fitting of the K —band image suggests an effective 
radius of ~ 11.5kpc. 

• A timescale of events can be built up from these results: star 
formation is triggered by some event; in this case, any tidal interac- 
tion or merger must have been relatively minor. AGN/radio source 
activity follows the onset of star formation by several million years. 

As a final consideration, we speculate on scenarios where the 
lack of evidence for a recent merger may be significant. As noted 
in the introduction, it is thought that the crucial ingredient for a 
radio-loud AGN is that the central black hole is rapidly spinning, 
a process which is most naturally explained by the coalescence of 
two smaller black holes (Wilson & Colbert 1995). Indeed, the cos- 
mic evolution of AGN activity can be well-described in terms of 
the expected halo merger rate (Mahmood, Devriendt & Silk 2005). 
However, the timescales involved in the black hole coalescence af- 
ter the halo merger are extremely uncertain; while it is thought 
that SMBHs can be brought within ~ one parsec of each other 
on a timescale similar to the dynamical friction timescale (e.g. 
Milosavljvic & Merrit 2001), the mechanism for the final merger 
is unknown. 

One possibility is that the formation of an accretion disc 
around one or both black holes provides the final necessary dynam- 
ical friction (e.g. Armitage & Natarajan 2005). If we assume that 
the central black hole of 9C J 1503+4528 was spun up as part of the 
AGN triggering process, the apparent simultaneity of the trigger- 
ing of both the optical AGN activity (requiring an accretion disc) 
and the radio jets (requiring the SMBH to be spinning rapidly) for 
9C J1503+4528 would demand that the final SMBH coalescence 
must take place on a very short timescale indeed. If this were not 
the case, and SMBH coalescence occurred more slowly, one would 
expect a period of radio-quiet AGN activity to precede the forma- 
tion of the jets, which is ruled out in the case of this source by the 
relative sizes of the EELR and radio source. 

Given the lack of evidence for a recent major merger, we there- 
fore propose that the host galaxy of 9C J1503+4528 contained a 
rapidly spinning SMBH prior to the present period of AGN activ- 
ity. In this case, the SMBH would have been spun up a significant 
time in the past, most likely during one or more mergers of smaller 
gas-rich discs which, in hierarchical clustering models, are the pro- 
genitors of ellipticals (e.g. Springel et al 2005). This scenario has 
the advantage that radio-loud AGN may be "re-lit" by relatively 
small gas injections from time to time, which provides a natural 
route for the necessary mechanical energy feedback in the intra- 
cluster medium (Nipoti & Binney 2005; Best et al 2005). However, 
both 9C J1503+4528 and the well-studied cluster radio sources are 
observed at relatively recent cosmic time. The time delay between 
the initial creation of a binary SMBH pair and the onset of jet activ- 
ity may have been significantly shorter in the early universe, where 
gas-rich mergers were significantly more prevalent. Such a scenario 
could be tested by investigating any time lag between the evolu- 
tion of the radio luminosity function at high redshift (which is still 
poorly constrained, e.g. Willott et al. 2001) and the halo merger rate 
as predicted by hierarchical clustering models. 
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Table A2. Changes to the best-fit YSP model with the inclusion of a power 
law component. Columns 1 and 2 give the percentage of the power-law 
component at a wavelength of ~ 3600A, and its spectral index, a (defined 
as Fx <x A Q ). Column 3 gives the best-fit old stellar population age in Gyr. 
The best-fit young stellar population age (in Myr) and mass (as a percentage 
of the total stellar mass of the host galaxy) are given in columns 4 & 5 
respectively. The reduced x 2 f° r the best-fit model is listed in column 6. 
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APPENDIX A: EMISSION LINE PROPERTIES OF 
9CJ1503+4528 

Here we present the results of our emission-line fitting for the paral- 
lel, perpendicular and combined 2-d frames. Fluxes relative to both 
[Oll]3727A and H/3 are tabulated in Table Al, as well as their 
equivalent widths. Although flux ratios quoted relative to H(3 are 
more common, the presence of Balmer absorption lines reduces the 
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Table Al. Spectroscopic properties of 9C1503+4528, determined from various spectra obtained perpendicular and parallel to the radio source. The spectra 
were extracted with various width apertures, so as to probe the properties of the central and exterior regions of the extended emission structures. The total 
integrated fluxes for the [Oil] and H/3 emission lines are provided, given in units of 10 _17 W m~ 2 . All other emission lines are quoted in terms of their flux 
ratios relative to both [Oil] and H/3 within the extracted one-dimensional spectra. This has been done due to the fact that the levels of Balmer absorption 
vary with position throughout the galaxy, and therefore a comparison with the H/3 line flux, though popular in the literature, may not be the most appropriate 
option for this source. The errors on the H/3 and [Oil] line fluxes are dominated by calibration errors, estimated to be < 10%. Errors on the other lines and 
this calibration error are added in quadrature with the errors arising from photon statistics. All quoted values in this table are corrected for galactic extinction 
using the E(B — V) for the Milky Way (Schlegel, Finkbeiner & Davis 1998) and the parametrized galactic extinction law of Howarth (1983). 



Details of extracted spectrum 


Combined, 2" 


Parallel, 2" 


Perpendicular, 2" 


Parallel, 5" 


Perpendicular, 5' 




[On] flux 


191.5 


227.5 


158.6 


248.0 


180.2 




H/3 flux 


11.7 


13.4 


10.0 


14.4 


10.2 


[Oll]3727 


[On] Flux Ratio 


100 


100 


100 


100 


100 




H/3 Flux Ratio 


1635 ± 177 


1701 ± 191 


1588 ± 195 


1715 ±231 


1766 ± 290 




Equiv. width 


182A 


192A 


170A 


150A 


123A 


[Nelll] 3869 


[Oil] Flux Ratio 


8.5 ± 0.9 


8.7 ± 1.0 


8.1 ± 1.0 


9.5 ± 1.1 


11.1 ± 1.5 




H/3 Flux Ratio 


139 ± 16.2 


148 ± 18.2 


128 ± 17.9 


162 ± 24.3 


196 ±36.7 




Equiv. width 


17A 


19A 


15A 


16A 


15A 


HC 3889 


[Oil] Flux Ratio 


0.7 ±0.2 


1.6 ±0.3 


1.1 ±0.3 


1.3 ±0.5 


< 0.3 




H/3 Flux Ratio 


12.1 ±3.4 


26.4 ±5.8 


16.7 ±5.6 


21.6 ± 7.9 


< 5.2 




Equiv. width 
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3.lA 


1.9A 
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0.6 ±0.2 


1.0 ±0.2 


< 0.3 
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2.7 ±0.4 


3.0 ±0.4 


2.7 ±0.5 


1.9 ±0.4 


2.0 ±0.7 




H/3 Flux Ratio 


43.7 ±6.5 


50.0 ± 7.3 


43.3 ±9.1 


31.7 ±7.0 


35.0 ± 13.5 
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3.9A 


2.3A 


2.lA 
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[Oil] Flux Ratio 


0.7 ±0.3 


0.9 ±0.4 


0.6 ±0.4 


0.8 ±0.7 


< 0.5 




H/3 Flux Ratio 


11.6 ±4.8 


14.4 ±6.3 


9.4 ±6.0 


13.6 ± 12.4 
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0.9A 


1.1A 
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[Oil] Flux Ratio 


3.3 ±0.5 


3.7 ±0.5 


2.8 ±0.5 


3.7 ±0.6 


1.5 ±0.5 




H/3 Flux Ratio 


54.2 ± 7.7 


63.0 ±8.7 
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5.7 ±0.9 
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100 


100 


100 


100 
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6.0A 
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10.7 ±1.1 
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23.7 ±2.5 
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30.5 ±3.2 
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383 ± 50.7 


523 ± 72.9 


446 ±83.1 
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38A 
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41.8A 
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usefulness of this ratio. The values in Table Al have been corrected 
for galactic extinction, using values for the Hi column density of the 
Milky Way taken from the NASA Extragalactic Database (NED), 
and the parametrized galactic extinction law of Howarth (1983). 



APPENDIX B: BEST-FIT YSP PARAMETERS IN THE 
PRESENCE OF A POWER-LAW. 

As part of our modelling, we consider the effects of the inclusion of 
a power-law component on the best-fit YSP parameters. The pre- 
ferred age of the young stellar population remains unchanged at 
5 Myr except for the strongest power-law contributions (> 25%). 
However, the reduced \ 2 values indicate that models with a strong 
power law contribution are not favoured. 



